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We study the effect of surface gravity waves on the motion of inertial particles in an 
incompressible fluid. Using the multiple-scale technique, we perform an analytical cal- 
culation which allows us to predict the dynamics of such particles; results are shown 
for both the infinite- and finite-depth regimes. Numerical simulations based on the ve- 
locity field resulting from the second-order Stokes theory for the surface elevation have 
been performed, and an excellent agreement with the analytical predictions is observed. 
Such an agreement seems to hold even beyond the formal applicability of the theory. We 
find that the presence of inertia leads to a non-ncgiigible correction to the well-known 
horizontal Stokes drift; moreover, we find that the vertical velocity is also affected by a 
drift. The latter result may have some relevant consequences on the rate of sedimentation 
of particles of finite size. We underline that such a drift would also be observed in the 
(hypothetical) absence of the gravitational force. 



1. Introduction 

The study of the Sto kes drift is a problem o f paramount importance both from a fun- 
damental point of view (|Van den Br oecl^'1999') and in connection with applicati ons, espe- 



cially in the area of sediment transport (Long uct-Higgins 1953; Nielsen 1992; Vitt ori &: Blondeaux 
1996t iBlondeaux. Brocchini fc Vittor i hooi iBlondeaux et aLll2012l ). As far as the first 



point is concerned, the Stokes drift is for instance responsible of important fiuid-mixing 
mechanisms such as mass and momentum transport near the free-surface, as well as ver- 



tical mixing enhancement owing to turbulent kinetic-energy production (jKantha et al 



In the ocean, the Stoke s drift is thought to be one im portant ingredient responsible 
for the Langmuir circulation (|Mc Williams. Sullivan fc Moe ng 1997). In relation to appli- 
cations, it is known that an accurate evaluation of the Stokes drift is important for the 
correct representation of surface physics in ocean general circulation models and ocean 
models at smaher sca les. Other relevant effects on the ocean circulation are discussed, 
e.g.. bv lNieisCTl (|l992[ ). 

Since the seminal paper by IStokes ( 1847 ). Stokes drift has been recognized as an 
important example that illustrates the d ifference between the Eulerian and the La- 
grangian statistics ( Longuet-Higgind 1986 ). It predicts that a fiuid particle (i.e. a tracer 
of negligible inertia) experiences a mean drift in the direction of wave propagation pro- 
portional to C/^/c, where U is the amplitude of the wave-induced velocity and c is 
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the wave phase velocity. Because the Stokes drift originates from the difference be- 
tween averages, it is relevant for all floating and suspended particles present in the 
water column, and not only for fluid particles considered in the original derivation. In- 
ertia of finite-size particles with density different from the fluid modifies Lagrangian 
averages with respect to Eulerian ones. This has important consequences on parti 
cle dispersion in both laminar and turbulent flows (see, e.g.. Squires fc EatonI 



1991 



Wang et a/."l992': 'Palkovich & Pumir"2004': Boffetta, De Lillo & Gamba"2004': 'Bee et al 



2006;,Martins AfonsQ2008: Martins Afonso, Mazzino k 011a 2009; Bee et a/. 2010 : Martins Afonso k Mazzino 



201C; Martins Afonso, Mazzino fc Muratore-Ginanneschi 2012i) . and we expect that in- 



ertia might affect the Stokes drift experienced by inertial particles. Previous studies in 
the field have investigated the case of parti cles close to be neutrah y buoyant in a velocity 
field generated by internal gravity waves (jGrinshpun et a/.ll2000l ) and small particl es of 
generic density in deep water in the presence of surface gravity waves (|Eamesll2008l ). 

Our main aim here is to push forward the analyses performed by these previous studies 
and to investigate the role of inertia on the resulting Stokes drift for linear and second- 
order Stokes waves over arbitrary depth. For this purpose, a multiple-scale expansion 
has been exploited to re-sum secular terms, the cumulative effect of which gives rise to 
a generalized Stokes drift. As a result of our analysis, we show that inertia induces a 
correction to the horizontal Stokes drift which is second order in particle inertia, and 
generates a vertical drift (a first-order effect) which modifies the sedimentation velocity. 
Interestingly, this vertical drift has a dynamical origin as it is active even in the absence 
of gravity, a remarkable result not pointed out in previous studies. The analytical results 
carried out by means of asymptotic methods are corroborated by a set of direct numerical 
simulations which extend the range of validity of our results beyond the perturbative 
regime. 

The remaining of this paper is organized as follow. In Section[2]we report the analytical 
calculations and results based on the multiple-scale expansion. Section [3] is devoted to 
the numerical results and their comparison with analytical predictions. In Section [4] we 
extend our results to the case of finite depth, and Section [5] is devoted to conclusions. 



2. Multiple-scale analysis of inertial-particle motion 

We consider the motion of small inertial particles transported by Stokes waves propa- 
gating in arbitrary depth h. The reason for considering a second-order expansion in the 
wave motion (Stokes wave) comes from the fact that the Stokes drift is a second-order 
effect in the wave amplitude. The two-dimensional, ir rotational and incompressible ve- 
locity field u = (u,w) induced by a Stokes wave is ( Whitham 1974t Kundu fc Cohen 
199(1 ) 

U cosh\k(z + h)] , 3U^ cosh[2k{z + h)] 

u = hr-n — - cos fca; - ujt) + '-^ cosMkx - ujt)] 

sinHkh) ^ ' 4c sinh^(fc/i) ^ ^ , ^ 

C/sinh[fc(z + /i)] . , W"^ smh[2k{z + h)] . 

w = . , , smikx - u)t) + — ^—^ '-^ sin 2 /ex - wt , 

sinh(fc/t) ^ ' 4c sinh4(fc/i) ^ ^ 

where x and z are the horizontal and the vertical coordinate, respectively; k is the 
wave number and uj is the angular frequency related to k via the dispersion relation, 
w = ^ gktwihikK) (strictly speaking, one should include the nonlinear correction to the 
dispersion relation; however, this turns out to be inessential in our analysis). The phase 
velocity is c = w/fc, and the maximum velocity \J at the surface (z — 0) of the first-order 
solution is related to the wave amplitude A by C/ = ujA. The equations ()2.1|) are obtained 
under the hypothesis of small steepness e = fcA, which, because of the relation between U 
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and A, is equivalent to the Fronde nnmber defined as Fr — U jc. Note that, in equations 
(|2.ip . the mean flow that should appear at the same order as the second harmonic is not 
included because we are dealing with a mono chromatic wave , and not with wave packets 



characterized by modulation length (see e.g. IWhithamlll974 p. 474 for a discussion) 



The motion of a small inertial particle transported by the fluid flow u through the 
Stokes drag (with typical response time r) and subject to gravity acceleration g is given 

by 



da; 

AV u-V ,^ Au 



V (2.2) 



+ (l-/?)g + /3— , (2.3) 



At T y i-J^ 



where x{t) and V{t) represent the particle position and velocity. In (|2.3p the added- mass 
effect has been taken into account via the dim ensionless number 13 — Zpi/ (pi-\-2p-o) , built 
from the fluid, pt, and particle, pp, densities ( Maxev fc Rilevl[l983t iGatignol 1983[) . 



In order to have an explicit expressi on for the parti cle velocity, we expand (j2.3p per- 



turbatively in r (fBalkovskv. Falkovich fc Fouxonll200ir) . to obtain: 



V = u + r{l-P){g-^^+r\l- /3)^ + 0{r') . (2.4) 

For simplicity, in the following we will first consider the deep-water limit, i.e. kh oo, 
for which the coefficients of the second order term in the velocity field vanish and (|2.ip 
simplifies to: 

u{x, z, t) = (Ue''^ cos{kx - ut), C/e*'^ sin{kx - wt)) . (2.5) 

In the case of deep-water gravity waves, the dispersion relation reduces to u — \fgk. 

Besides the steepness, another dimensionless number, F = (1 — /3)gTc/U'^, can be 
introduced as the ratio of the bare gravitational sedimentation velocity, (1 — (3)gT, over 
the drift velocity, U^/c. Note that this number can be written as F = (1 — /3)St/Fr^, 
where we have introduced the Stokes number St = ujt. Indeed, since the velocity field 
(|2.5|) decays exponentially with z, the effect of the Stokes drift is modulated by the 
sedimentation process. Therefore, in the perturbation expansion, one has to decide at 
which order in e the gravitational term has to be taken into account. We have analyzed 
different, physically relevant, cases and we found that the most interesting regime is when 
F = 0(1). Physically, this choice corresponds to a situation in which gravity is relatively 
small and particles experience a finite horizontal Stokes drift during sedimentation. In 
other words, the particle settling velocity in still fluids, (1 — l3)gT, made dimensionless 
via the wave velocity c, is assumed to be asymptotically cx U^/c^, i.e. O(e^). 

Upon momentarily introducing the dimensionless variables x kx , t i—i' ujt and u 
u/U, we can expand u and its Lagrangian derivatives at the second order in e, and 
substituting in (|2.2p - (|2.4p we obtain for the particle motion: 

x = e[u~ St(l - P)w - St2(l - /3)ii] + e2St^(l - I3)e^^ + 0{SV^) + 0{e^) , (2.6) 

z = e[w + St(l - /3)m - St^(l - /3)w] - [F + St(l - ^)e^^] + 0{St^) + 0{e^) . (2.7) 

A standard expansion of the coordinate as x — Xq + exi + e^X2 + ... leads to sec- 
ular terms in the equation at the second order. Therefore, we are led to analyze the 
problem within the multiple-scale formalism. To this aim, we introduce an independent 
slow time variable T = e^t — the specific choice of the power 2 as the exponent of 
e is obviously dictated by the fact that the secularity appears in the equation at the 
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second order of the expansion in e — , and we consider the dyna mical variables as a func- 
tion of both times, x(t,T) = xo(t,T) + exi(t,T) + e'^X2(t,T) (|Bender fc Orszadll978 



Bensoussan. Lions fc PapanicolaoTjll978l ). At the second order, from the solvability con- 
dition which removes secular terms, we obtain an expression for the mean Lagrangian 
velocity which, going back to dimensional variables, gives the drift velocities: 

Ud = — \l~ /?(! - /3)Se] e2M-o-(i-«9rt] ^ (2.8) 
c 

= -(1 - (3)gr - (1 - /3)St— e2'=["«-(i-'5)frt] . (2.9) 

c 

In the limit of fluid particles, St = and /3 = 1, the above expressions recover the drift 
derived by Stokes: Ud = e^'^^''t/^/c, Wd — 0. Inertia induces a correction of order St^ 
to the horizontal drift velocity (I2.8p . The interesting result is that inertia produces a 
drift velocity also in the vertical direction, which corrects the sedimentation velocity. 
The vertical drift is present even in the absence of gravity (5 = 0), and is a consequence 
of the interplay between delayed particle dynamics (due to inertia) and asymmetric 
vertical dependence of the velocity field. The vertical drift (|2.9p can be rewritten as 
Wd = —(1 — /3)5t{1 + Fr^e^'^'I^"^^^"^)^'^*!}, from which we see that the correction to the 
bare sedimentation due to gravity is simply proportional to Fr^. To our knowledge, the 
correction to sedimentation velocity induced by water waves has never been discussed 
before. We also notice that, here and in what follows, no long-time divergence appears for 
light particles (/3 > 1, i.e. F < 0), because in that case the particle evolution is bounded 
by the water surface (z = plus wave corrections). 

Another simple, and physically relevant, prediction one can derive from (|2.8p is the net 
displacement of particles from the point of release. For heavy particles (fi < 1) released 
at the surface {zq = 0), a time integration of (I2.8P from to 00 gives for the total 
displacement: 

2fc(l - /3).grc L ^ 2k (1 - /3)St ^ ' 

Of course, expression (|2.10p is in principle valid only for small St, such that Aa; > 0. 



3. Numerical simulations 

In this Section we report the numerical results obtained from the integration of (|2.2p 
and (12.31) . with the aim of verifying the analytical predictions of the previous Section, 
and also to check the robustness of these predictions with respect to the expansion 
parameters. 

Figure [T] shows two typical examples of trajectories of slightly heavy and light particles 
in the deep-water limit. From the trajectories the drift velocity is obtained by computing 
averages over the Lagrangian period (which is different from the Eulerian, wave period). 
An example of the resulting velocity components is shown in Fig. [2] for a wave with 
e = 0.33. Despite the fact that the steepness is not so small, the agreement with the 
theoretical prediction (|2.8p and (|2.9p is very good. From the plot of the vertical velocity 
Wd we see that the relative correction induced by waves to the bare sedimentation velocity 
at the initial time is Fr^ « 0.11. 

For steeper waves, the agreement between numerical simulations and perturbative 
predictions worsens, but (|2.8p and p.9p still capture qualitatively the fact that bigger 
waves induce larger corrections to the sedimentation velocity. One should also realize 
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Figure 1. Two examples of trajectories of slightly heavy (/3 = 0.99, continuous line) and light 
{(3 = 1.01, dotted line) particles transported by a linear wave 1)2. 5|) with e = Fr = 0.33 and 
St = 2.22 in deep water. The initial position for particles is a;(0) = 0, 2:(0) — (heavy) and 
x(0) = 0.13, z{0) = -0.5 (light). 



that for a larger steepness other effects such as wave breaking, clearly not included in 
the theory or simulation, can take place. 

The total horizontal displacement of heavy particles released at the surface of deep 
water for different values of the parameters is shown in Fig. [3l together with the predic- 
tions given by (|2.10p . The agreement is very good not only in the perturbative regime of 
small St in which (|2.10l) is derived. Deviations are observable only for St = 0(1). 

As mentioned in Section [21 the analytical result (|2.9I) implies that inertial particles 
have a vertical drift even in the absence of gravity. This drift acts in the same direction 
of sedimentation as it has the same sign of the gravitational term, independently on /3. 
By comparing (|2.8p and (|2.9p . we see that for g = the mean motion is along a straight 
hne, with slope given by w^/ud = —(1 — /3)St/[l — — /3)St^]. A simple consequence 
is that in this case the total displacement of a heavy particle over deep water diverges, 
in agreement with ()2.10|) with 5 — > 0. Figure |4] shows the slope of the mean motion as 
a function of St for different values of (3. Again, for small and moderate values of the 
parameter, the agreement with the analytical prediction is very good. 



4. Predictions for arbitrary depth 

In this Section we generalize the calculations of Section [2] to the velocity field (|2.ip 
induced by a Stokes wave in water of arbitrary depth h. The calculation follows the same 
steps of the deep-water limit, with the additional complication due to the presence of 




Figure 2. Numerical (squares) vs. theoretieal (solid line) drift velocity: horizontal (upper panel) 
and vertical (lower panel) components. The dashed line represents the settling velocity in the 
absence of wave motion wa = —(1 — P)gT. Parameters: e = Pr = 0.33, St = 2.22, /3 = 0.99. 
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Figure 3. Total horizontal displacement of heavy particles with /? = 0.9 (squares), /3 — 0.99 
(circles) and P = 0.995 (triangles), settling beneath a linear wave with e = Fr = 0.33, as a 
function of St. Lines represent theoretical predictions (|2.10|) . 



new terms. The final result, which generalizes equations (|2.8l) - (|2.9p . reads 

C/2 1 -/3(l-/?)S t^ 
2sinh^(fc/i) 



= „ . 2 cosh{2fc[zo + h-{l- p)gTt]} , (4.1) 



= -(1 - l3)gT ^ . j; sinh{2fc[zo + h - [1 - /3)gTt]} . (4.2) 

c 2sinh (fc/i) 

We first notice that, by taking the limit kh oo, the above expressions recover the deep- 
water result of Section O We also observe that, remarkably, the second-order term in the 
Stokes wave expansion (i.e. the one proportional to in (12.11) ) does not contribute in 
the drift (|4.ip - (l4.2p . as long as we stop the expansion (I2.6p - (l2.7p at the second order in 
e. Its effect will eventually appear in higher-order terms in the expansion, which are not 
considered here. From this point of view, the perturbative analysis for finite depth has 
the same structure of the infinite-depth limit. 

We have performed numerical simulations of particle motion given by (|2.2|) and p.3p 
in the velocity field (|2.1|) . Also in this case we have found a very good agreement between 
the numerically-evaluated drift velocities and the analytical predictions (|4.ip and (|4.2p . 
even for values of the expansion parameters relatively large, confirming the robustness 
of the multiple-scale technique. 



5. Conclusions 

We have considered the problem of Stokes drift induced by water waves on small 
inertial particles with two complementary perspectives. On the one hand, our results 
give the correction to the horizontal Stokes drift induced by inertia. This correction is 
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St 

Figure 4. Ratio of vertical to horizontal drift velocity Wd/{{1 — P)ud as a function of St 
for heavy particles in the absence of gravity {g = 0). Parameters: e = 0.33, /3 = 0.9 
(squares), /? = 0.99 (circles), /3 — 0.995 (triangles). The Une represents the prediction 
Wd/{1 - fi)ua = -St/[1 - P{1 - P)St\ 

found to be second order in the particle Stokes number, with a sign which depends on 
the particle density relative to water. On the other hand, we also obtain a vertical drift, 
which therefore represents a correction to the sedimentation velocity induced by wave 
motion on the surface. This effect, which results to be first order in the Stokes number, 
has never been discussed before and is of possible relevance, e.g., in the field of sediment 
transport in coastal regions. 

We conclude observing that, although the present analysis is performed in the ideal 
world of linear two-dimensional water waves and in the absence of any interactions with 
physical boundaries, we expect that our main findings will survive in more complex and 
realistic situations. Indeed, we have already shown that, in the more complex case of 
second-order waves on arbitrary depth, our results do not change qualitatively. There- 
fore, it would be extremely interesting to study the drift of inertial particles, and its 
effect on sedimentation, in more realistic simulations of wave motion and in laboratory 
experiments, where a precise determination of the mean velocity and falling velocity is 
possible. 
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